Geometrical Frustration and Static Correlations in a Simple Glass Former 
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We study the geometrical frustration scenario of glass formation for simple hard sphere models. 
We find that the dual picture in terms of defects brings little insight and no theoretical simplifi- 
cation for the understanding of the slowing down of relaxation, because of the strong frustration 
characterizing these systems. The possibility of a growing static length is furthermore found to be 
physically irrelevant in the regime that is accessible to computer simulations. 
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The ubiquitous glass formation and jamming still puz- 
zle physicists. How can molecular and colloidal sys- 
tems slow down so abruptly without obvious structural 
changes? In response to this conundrum, theoretical ap- 
proaches inspired by spin-glass physics have long postu- 
lated a role for a "hidden" static length associated with 
the dynamical slowdown. Following the rigorous iden- 
tification of a growing point-to-set length accompanying 
diverging relaxation times in structural glasses [l|- 1| , a se- 
ries of "order agnostic" proposals for static correlations 
in supercooled liquids have also recently flourished 
and their analysis is ongoing. Yet specifying a relevant 
amorphous order parameter that captures these mate- 
rials' rich phenomenology while providing geometric in- 
sights into the underlying microscopic mechanism is still 
sought after. A scenario for growing geometrical order 
proposed some time ago by Sadoc and Mosseri [3] as well 
as by Nelson and coworkers [|| is often considered by 
many to suit this purpose. In addition to encouraging 
the enumeration of preferred local structures, e.g., [9f, it 
has indeed led to the development of a theoretical appa- 
ratus for the glass transition based on geometrical frus- 
tration [1, E3]. Yet, in spite of its marked impact on the 
structural analysis of dense fluids, this proposal remains 
largely untested in three-dimensional (3D) systems. 

Geometrical frustration is canonically illustrated by 
considering the behavior of spherical particles of diameter 
a. Because regular simplices (triangles in 2D, tetrahedra 
in 3D, etc.) are the densest possible local packings of 
spheres, they are expected to play a central role in liquid 
organization, e.g., [ill EH- I n 2D Euclidean space, in- 
teresting physics results from the fact that simplices can 
assemble into the triangular lattice |8j, and spatial cur- 
vature frustrates the regular assembly of disks (l3j ]. For 
Euclidean space in dimension d > 3 simplices cannot tile 
space without defects, but in d = 3 they can form per - 



results in irreducible disclinations that sit on disk cen- 
ters and for which the coordination obtained by a De- 
launay decomposition differs from six. Similarly, in 3D 
each edge between nearest-neighbor pairs is shared by five 
other tetrahedra; flattening space generates disclination 
lines of "bond spindles" that are shared by q ^ 5 tetra- 
hedra. Periodic arrangements of these disclinations form 
the complex crystal structures known as Frank-Kasper 



phases 15j. Yet even in amorphous configurations at 



feet icosahedra on a relatively gently curved space (14 1. 
The defects that result from uncurving this singular space 
back to the Euclidean variety can be understood by di- 
mensional analogy. Each particle in a perfect 2D trian- 
gular tiling of disks is part of six triangles. Curvature 



small frustration, the simple Voronoi polyhedra that ac- 
commodate the presence of spindles with 9 / 5 pro- 
vide topological constraints for the propagation of defects 
from one particle to the next, which results in disclina- 
tion lines [8J. A denser fluid, in which the proportion of 
q = 5 spindles grows and conversely that of q ^ 5 spin- 
dles shrinks, should thus see disclination defects play a 
more important role. Disclination lines passing one an- 
other correspond to activated events, possibly affected by 
topological constraints Q ; the theoretical framework sug- 
gests a causality between the dynamical slowdown and a 
growing static, structural correlation length underlying 
the fragility of the glass-forming fluid. In this letter, we 
critically examine this proposal and find that, in spite 
of its elegance, it does not hold for the system to which 
it is more directly related, i.e., simple 3D hard spheres. 
Through a variety of measures of static order, we also 
consider alternate definitions of correlation lengths and 
explore in what regime a growing static order could rea- 
sonably be associated with a dynamical slowdown in the 
regime accessible to computer simulations. 

Testing this intrinsically geometrical theory on iden- 
tical 3D hard spheres is problematic because nucleation 
interferes with the slowdown, so the crystallization drive 
must be reduced. In case the fluid structure were 
to sensitively depend on the nature of these perturba- 
tions, we do so in two different ways: (i) a 50:50 hard 
sphere binary mixture with a 1.4:1 diameter ratio whose 
glass-forming properties have been extensively character- 
ized Til- 18 1 . and (ii) a mixture of hard spheres with 
the smallest non-crystallizing diameter polydispersity, 
8.5% 19]. The average number of tetrahedra wrapped 
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FIG. 1: Evolution of q and (left inset) of the probability dis- 
tribution p(q) with density. The statistical honeycomb limit 
is given for reference, (right inset) Network of q = 6 spindles 
(rods) at <j> = 0.58. 



FIG. 2: Long-time limit of the overlap, (inset) Time t evo- 
lution of the overlap for c p i n = 0.11, 0.45, 0.68, 0.74, and 
0.79, from left to right, in (/3mcr 2 ) -1 / 2 units with mass m, 
the larger particle a and inverse temperature /3 set to unity. 



around a bond, q, lies within two simple limits. First, 
all finite-density configurations should have fewer tetra- 
hedra per spindle than a Poisson process (an ideal gas), 
where q = 144tt 2 / '{24ir 2 + 35) w 5.228 Q. Second, al- 
though in curved space the optimal number of simplices 
per bond can be as low as q = 5, in Euclidean space a 
more stringent limit q = 2ir/ arccos(l/3) s» 5.104 is ob- 
tained from the fictitious "statistical honeycomb" con- 
Figure Q] shows that both densifying 
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struction 

fluids present a growing polytetrahedral character. The 
average spindle coordination decreases, seemingly toward 
its optimal value, with increasing packing fraction <f> for 
both models, as does its distribution p(q) (Fig. Q}. The 
growing simplex order is also quite different from that ob- 
served in the face-centered cubic crystal phase. Because 
the structural properties of the two models are robustly 
similar, we only consider (i) for the rest of the analysis. 

Surprisingly, even for the densest systems equilibrated 
the disclination network remains highly branched, with 
multiple defect lines stemming from each vertex. The 
inset of Fig. Q] illustrates this situation for the q = 6 
spindle network. The typical spacing between defect 
spindles £dcfcct = c d e fect> usm § a defect concentration 
Cdefect = Ylq c q(l ~ 5) 2 that puts more weight on higher- 
order defects, indeed grows by no more than 1 — 2% over 
a density range over which the relaxation time goes up 
by several decades. Extrapolating the results to higher 
densities using the statistical honeycomb limit further in- 
dicates that the growth of ^defect remains small over the 
entire accessible amorphous regime <j> < 0.65. Actually, 
polytetrahedral order is bound to saturate as a result of 
the intrinsic frustration of Euclidean space. The satura- 
tion length corresponding to the maximal spatial exten- 
sion of simplex order estimated from the radius of sphere 
inscribing the gently curved space from which this argu- 
ment derives is ~ 1.59ct 14[. Alternatively, the typical 



distance between defects in an ideal tetrahedral structure 
threaded only by q = 6 spindles is only ^defect — 0.99a, 
although it is worth noting that the average distance be- 
tween spindles itself is but ~ 0.3c. A similar result is 
obtained for the spatial correlations associated with frus- 
trated local order through an analysis of the spatial de- 
cay of the bond-orientational order correlation function 
Ge(r) [22|, Hj|. No matter how it is precisely defined, 
the associated correlation length £ 6 does not increase by 
more than a few percents. Even in 3D, hard spheres are 
therefore sufficiently frustrated to make the dual picture 
of amorphous particle packings in terms of spindle defects 
rather uneconomical at all densities. 

In order to remove any possible doubt as to whether al- 
ternative static lengths due to tetrahedral or other order 
types are present or not, we turn to the order-agnostic 
penetration length £ p 2J, |25j , which, like the point-to- 



set length £ps [4( , characterizes the influence of boundary 
conditions and is expected to diverge with the relaxation 
time 0]. It is obtained by pinning a random selection 
of particles from an equilibrated fluid configuration, and 
measuring the overlap between the initial and final con- 
figurations after a long time t has elapsed, 



Q(c) = lim 

t— >OC 



E>,(QK(<)) 



(i) 



where m is the occupancy of a spatial cubic cell whose 
volume is similar to that of the smaller particles in the 
system in order to prevent multiple occupancy [24j . Sub- 
tracting the random overlap contribution (m) leaves a 
quantity that grows from low to high as the pinning 
concentration c P i n increases (Fig. |5J); the crossover is 

£ p ~ Cpjj/ 3 . Operationally, we define £ p as the value 
of the average distance for which the overlap falls below 
0.18 (Fig. |2J| . The extracted length is not very sensitive to 
this choice, provided it is intermediate between low and 
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FIG. 3: Various static lengths rescaled to unity at <p = 0.40 
(£o) together with the lower bound from Eq. ([2]) Abound and 
the dynamic length [It}]. (inset) Measures of dynamical slow- 
down. Lines are guide for the eyes. 



high overlap. We stress that focusing solely on the low- 
overlap regime provides no information on £ p as it only 
depends on the standard pair correlation function and 
therefore on trivial two-point correlation lengths £2 • This 
result, which can be checked explicitly by considering the 
linear response to a vanishingly small c p - ln , remains true 
so long as one remains in the low overlap perturbative 
regime 



261 ] . It also casts some doubt on the relevance of 



a recently proposed scaling 27] . where the observed lin- 
ear dependence on concentration suggests instead that 
only trivial static lengths are probed. 

Figure [3] shows that the penetration length £ p increases 
only very modestly over the dynamically accessible den- 
sity range. For sake of comparison, we display in Fig. [3] 
an estimate of £2 evaluated from the pair correlation 
function g(r), which, as is well known, changes only 
slightly with density. The penetration length increases 
only about 20% more than £2- (Artificially breaking 
down structural order into radial and orientational con- 
tributions suggests that the latter is at most comparable 
to the former over the density range studied, further sup- 
porting the spindle analysis.) Note that in view of the 
small variation of all the static lengths, ^defect, ^6, £2 
and £ p , trying to devise a crisper measuring procedure 
is unnecessary as it will not qualitatively alter the con- 
clusions. Strikingly, the "dynamic length" £dyn charac- 
terizing the spatial correlations in the dynamics and as- 
sociated with dynamical heterogeneities grows markedly 
over the same density range. Whereas the change in the 
static length is measured in fractions of their low-density 
value, £rivn grows by a factor of almost 7 when reaching 
4> = 0.59 [lTj. The diffusivity D and the structural re- 
laxation time t meanwhile change by about 4 orders of 
magnitude (Fig. [3]). 

Although the above results may come as no surprise to 
those who believe the dynamical slowdown to be a purely 
kinetic phenomenon involving no growing static length 



scale, it is nonetheless worth checking whether one does 
not violate the bound between relaxation time and static 
correlation length put forward by Montanari and Semer- 
jian [H , t < To exp (i? £p S ) , where To is a constant setting 
the microscopic time scale. The coefficient B depends on 
density (or temperature for a glass-forming liquid) and 
is such that when £ps ~ a the right-hand side describes 
the "noncooperative dynamics" of the model 0]. Using 
an Arrhenius-like argument for activation volumes [28j . 
we note that, all else being equal, higher pressures triv- 
ially rescale the free-energy landscape and thereby slow 
the dynamics. For a hard-sphere fluid, one then expects 
B cx f3P where P is the pressure. It should be stressed 
that the upper bound of r diverges with the pressure even 
in the absence of any growing £pg , as when approaching 
T = for an Arrhenius temperature dependence. In the 
low and moderate density fluid, the relaxation time in- 
deed follows r((j)) ~ T\ ow ((j)) — t exp{K(3P((j))] with K a 
density- independent constant. One then finds that 



log[r(0)/r o ] 
log[ri ow (0)/r o ] 



1/3 



< 



£ps,o 



(2) 



where £ps,o is the low-density limit of £ps- Equation ^ 
thus provides a lower bound for the growth of a static 
length imposed by the dynamical slowdown. 

To assess whether the above bound is satisfied or not, 
one needs an estimate of £ps- The direct approach would 
be to consider the effect of pinning the boundary of a 
spherical cavity on the fluid inside, but one may reason- 
ably expect that the penetration length studied above 
gives a rough estimate of £ps ■ Near a random-first-order 
transition or near any first-order transition £ps ~ £p , 
but far from such transitions, which is the case studied 
here, one expects £ps ~ £ P - In any case, as seen in Fig. [3J 
the bound given by Eq. ([2J increases only slowly in the 
dynamically accessible domain and is already satisfied by 
£ p . Note that this moderate growth of the bound further 
illustrates that hard spheres are not in fact very "fragile" 
in the regime up to <j> — 0.59, showing only a limited de- 
viation with respect to the low-density behavior, which 
is in line with what is found for other simple fluids, such 
as the Lennard- Jones glass- forming liquids |3(J. These 
observations may well correlate with the fact that most 
3D fluids of spherical particles are strongly frustrated in 
the sense discussed above. 

These results indicate that the growth of a static length 
is not the controlling factor behind the relaxation slow- 
down in the range of density considered. This finding 
points to a mechanism for the slowdown that is either 
essentially "noncooperative", or akin to that predicted 
by the mode-coupling theory (MCT) [H [13]; in both 
cases, the growth of a dynamic length is not accompa- 
nied by that of a static length. We cannot, however, 
draw any general conclusion on this question beyond this 
regime. In thermodynamic- based theories [Iol[33| . it is at 



4 



these higher densities (or lower temperatures for a liquid) 
where cooperative behavior becomes dominant, and the 
dynamical slowdown is predominantly controlled by the 
growth of a static length. This regime is unfortunately 
mostly beyond present-day computer resources [34j. A 
modest indication that a crossover takes place may, how- 
ever, be given by the data for £ p and the bound, which 
both appear to display a steeper increase near <f> = 0.59. 

We have first shown that 3D hard spheres, like d > 3 
hard spheres and many 3D simple glass formers, are too 
strongly frustrated for the dual picture in terms of de- 
fects [8[ to bring any useful simplification. One may 
wonder if there exist other liquids with a different type 
of locally preferred order for which frustration is weaker 
and the picture can be put to work. Whereas this frustra- 
tion regime can be achieved in 2D by curving space [l3| , 
no such clear-cut example of simulation-accessible glass 
formers in 3D Euclidean space has yet been devised. Sec- 
ond, we have shown that within the regime accessed here 
static lengths grow very slowly, yet in a way that is com- 
patible with the bound recently put forward between re- 
laxation time and static length. For systems similar to 
hard spheres, these results severely constrain the type of 
ordering that can develop and place serious doubts on the 
pertinence of local-order analysis in the moderately vis- 
cous dynamical regime. The dynamic length's significant 
increase points instead to a decoupling between the in- 
creasingly heterogeneous character of the dynamics and 
its cooperative origin in terms of structural or thermo- 
dynamic quantities. A challenge would be to search for 
a possible crossover at still higher densities. 
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APPENDIX: GEOMETRICAL TILINGS OF 
SIMPLICES 

We review the regular polytopes and hyperbolic cover- 
ings that correspond to defectless tilings of simplices on 
spaces of constant curvature. In particular, we aim to 
show the following. 

1. A regular tiling of tetrahedra is found on the sphere 
that inscribes the remarkably large 4D platonic 
polytope {3,3,5}. The inscribing sphere's curva- 
ture, 7r 2 /25, is notably much smaller than that of 
the generalized octahedron {3,3,4}, 7r 2 /4. 

2. By contrast, on the d-dimensional sphere with d > 
3, {3 d_1 , 4} is the only way to obtain a regular tiling 
of simplices resulting in large geometrical frustra- 
tion of the corresponding Euclidean space, which 
hints at a possible limit role of simplex ordering in 
3D. 

3. The regular simplicial tiling {3,3,3,6} is found on 
a hyperbolic space of curvature larger (in absolute 
value) than the curvature of the sphere that in- 
scribes the tiling {3,3,3,4}, while the hyperbolic 
tiling {3,3, 6}'s infinite edge length disqualifies it 
as a possible reference. 

4. The homogeneous statistical honeycomb {3, 3,q} 
gives is an upper limit for the density of simplex 
tiling in d — 3, and a lower limit to q. 

Notations 

We first establish some notation. The Euclidean space 
of dimension d is denoted R d and the hyperbolic space of 
dimension d is denoted H d . The sphere (or hypersphere) 
of (surface) dimension d is the set of all vectors of length 
1 in R d+1 and is denoted S . Hence S 2 is the standard 
radius one sphere in R 3 we all know and love. 

Following Schlafli's notation (see [l|, Chapter 7]), 
the polytope {a\,...,ad} is composed of polytopes 
{ai, . . . , a,d_i} with vertex figure {02, . . . , aj}- The ver- 
tex figure of a d-dimensional polytope is a (d — 1)- 
dimensional polytope whose vertices are obtained by tak- 
ing the middle of each edge emanating from a vertex. 
This inductive definition deserves an illustration. The 
basic example is {p}, the regular polygon of p sides in 
R 2 . Next we have the regular polytopes in R 3 . Those 
are represented by the symbols {p, q} for a polytope hav- 
ing faces made of {p}s and vertex figures made of {q}s. 
Note that because the polytope is regular, one can in- 
scribe all the vertices on a sphere and thus {p, q} gives a 
tiling of the sphere S 2 . For instance, the cube {4,3} is 
composed of squares (the regular polygon {4}) and the 
vertex figure is {3}, a triangle. The cubic tiling {4,4} in 
R 2 is composed of 4 squares surrounding each vertex. 



Throughout the text, we use a measure q of the num- 
ber of tetrahedra sharing a particular edge between 
nearest-neighbor pairs and q for the average amongst 
all pairs of nearest-neighbor. This quantity generalizes 
to higher-dimensional setting in the following manner: 
in a (i-dimensional space (be it M. d ,S d or H d ), a set of 
(i + 1) points all nearest-neighbors to each other form a 
i-simplex (generalized tetrahedron). A 1-simplex is a line 
segment, a 2-simplex is a tetrahedron, etc. As in three- 
dimensional system, the dimension (d— 2)-simplices play 
an important role. We call them spindles. Given a par- 
ticular (d — 2)-simplex, we let q measure the number of 
d-simplexes that contain it. Obviously, q is the average 
of q over all spindles. 



Schlafli's notation and spindles 

Of primary importance for this paper is the rela- 
tion between Schlafli's notation and the number q. Let 
{ai, . . . , a,i} be a regular tiling in R d or S d . One can show 
that q = ad- Instead of exhibiting the proof of this result, 
we give an example. Consider the simple cubic tiling of 
R d , given by {4, 3 d " 2 , 4}. Here, a (d - 2)-dimensional 
hypercube is given, for example, by x\ — x-i = and 
< xi < 1 for 3 < i < d. The various d-dimensional hy- 
percubes containing this (d — 2)-dimensional hypercube 
are given by ±2:1, ±X2 € [0, 1], where the signs are chosen 
independently. Hence 9 = 4. 

Note that any polytope {a±, ...,a,d} in R d+1 can be 
seen as inscribed in a sphere S d whose center is the center 
of mass of the polytope. By inflating the polytope to 
make it round, we see that it produces a regular tiling 
of S d . For instance, the hypercube {4, 3 d-2 } in R d is 
also a tiling of S' d_1 . By looking at the usual case {4, 3}, 
one sees that it gives q — 3 on S 2 . The same argument 
illustrates that for the tiling {4, 3 d " 1 } of S , one also has 
g = 3. 



Various tilings in various spaces 

Consider the case of the generalized octahedron, 
{3 d " 1 , 4}, in R d+1 . It gives a regular tiling of S d with 
q = 4. This tiling is dual in R d+1 to {4,3 d " 1 } which 
can be thought as the cube [—1, l] d+1 . So the vertices of 
the dual generalized octahedron are ±e*j for the standard 
basis {ei, . . . , ed+i} of R d+1 . Obviously, the distance be- 
tween two vertices is \/2 via straight line, or ~ in S d . 
To obtain a unit distance in S , one has to rescale the 
sphere to be of radius — , hence the sphere has a curvature 
^ ~ 2.47. 

As pointed by d p.20], the 120 vertices of the 600- 
cell {3, 3, 5} all belong to the hypersphere S 3 with radius 
equal to the golden ratio (r = 1+ 2 V ^ ) if the edges are of 
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Tetrahedral honeycomb 


in 


curvature 


{3,6} 


R 2 





{3, q} with q > 7 


H 2 


— 4arccosh 2 (2 _1 esc -) 


600-cell {3, 3, 5} 


S 3 


tt 2 /25 ~ 0.39 


octahedron {3 n_1 ,4} 




tt 2 /4 ~ 2.47 


{3,3,3,6} 


H 4 


— 4arccosh 2 r ~ —4.505 



TABLE I: Summary of tiling results with the curvature re- 
ported setting (7 = 1. 



unit length. The alternative expression r = icsc(j^) is 
useful for our computation. Indeed, imagine a isosceles 
triangle whose equal sides measure r and opposite side 
measure 1. The angle opposite to the side of measure 1 
is 



9 = 2 arcsin 



2 arcsin 



in ( sin( — ) ] = — . 
V V 10V 5 



So, if one desires the spherical distance between vertices 
in the {3,3,5} regular tiling of S 3 to be 1, the radius 
must be — ^rs = — • The curvature of this sphere is thus 

T7T/5 7T " 

4 =s 0.39. 



The tiling {3,3,6} in hyperbolic 3-space H 3 is some- 
what of a degenerate case. Despite the cells being of 
finite volume, all of its vertices are at infinity hence the 
edge length is infinite, which makes the curvature im- 
possible to rescale; see 0, p. 202]. In if 4 however, the 
tiling {3, 3, 3, 5} is perfectly well-defined, with finite dis- 
tance between vertices. According to [3J, p. 204 and table 
p. 213 ], the edge-length is 2(f) and cosh</> = r, hence the 
edge-length is 2arccoshr ~ 2.12. The expressions ob- 
tained by Coxeter are consistent with the results in stan- 
dard hyperbolic space of curvature —1. Rescaling the 
space to make the edge-length 1 changes the curvature 
to — (2arccoshr) ~ —4.505. 

Still following 0, Chapter 10], one sees that the regu- 
lar tilings in H 5 are not tetrahedral and that there are 
no regular tilings in hyperbolic space of six or more di- 
mensions. In H 2 , any {3,g} with q > 6 is possible, 

with edge-length 2</> where cosh</> = C ?*J = ^- 3 -; d, 

p. 201]. Rescaling to get edge-length 1, one gets a curva- 
ture — 4arccosh 2 (2~ 1 esc £). 

The results for d > 3 are summarized in Table fl] Note 
that in the H 2 case, the packing of an increasing num- 
ber of hyperbolic triangles around a vertex, keeping the 
edge- length 1, makes the curvature go down from approx- 
imately -1.19 (for {3, 7}), to -2.34 (for {3, 8}), to -3.44 
(for {3,9}), with limit — oo as q — > oo. 



Statistical honeycomb 

Coxeter suggested that nature tries to approximate 
regular tilings {p, 3,3} for p between 5 and 6 in soap 
froth 0, Chap. 22]. Obviously, the soap froth itself pro- 
duces an irregular tiling, but one could look at the aver- 
age structure, the statistical honeycomb. Taking its dual, 
as suggested by [f| , we consider the corresponding statis- 
tical honeycomb {3, 3, q} produced by wrapping as many 
tetrahedra as possible around a given spindle, even if that 
means wrapping a fraction of a tetrahedron. Because the 
angle formed by two adjacent edges in a perfectly regu- 
lar tetrahedron is arccos(l/3), the ideal number is thus 

As proposed in Ref. [5j, we can consider an ideal con- 
figuration where every bond belongs to q tetrahedra. In 
the Voronoi decomposition corresponding to this config- 
uration, the cells have Z identical faces corresponding to 
the Z nearest neighbors of a particle, each of them being 
a polygon with q faces. This Voronoi cell can be seen as 
living on the surface of a sphere. If it has V vertices and 
E edges, we must thus have V — E + Z = 2, using Euler's 
relation. Counting the pairs (edge, face), with edge adja- 
cent to face, in two different ways, we get 2E = Zq, and 
counting the pairs (edge, vertex) in two different ways, 
we get 3V = Zq. Overall, we obtain that in this ideal 
structure, each particle has a number of neighbors 



Z = 



12 

6 — q 



13.4. 



Note that in this construction, all quantities are also by 
definition their own average. 

Now one can compute the volume occupied by that fic- 
titious Voronoi cell, and obtain, as did Ref. [HJ , a volume 
of 



Za 



area(q-sided regular polygon) 



3 2 

_ Za 3 gcot(|) 
144 

The packing fraction is then given by 

, ^sphere 



HI) 3 



coll 



Vc 



0.7796. 



cell 



It should be noted that this result is also Roger's bound 
on the maximal packing fraction for any sphere packing 
in 3D, not only simplex tilings. In all cases, it gives 
an upper bound on the density of simplex tilings in 3D 
systems. 
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